A quantum liquid with deconfined fractional excitations in three dimensions 
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Excitations which carry "fractional" quantum numbers are known to exist in one dimension in polyacetylene, 
and in two dimensions, in the fractional quantum Hall effect. Fractional excitations have also been invoked 
to explain the breakdown of the conventional theory of metals in a wide range of three-dimensional materials. 
However the existence of fractional excitations in three dimensions remains highly controversial. In this Letter 
we report direct numerical evidence for the existence of an extended quantum liquid phase supporting fractional 
excitations in a concrete, three-dimensional microscopic model — the quantum dimer model on a diamond 
lattice. We demonstrate explicitly that the energy cost of separating fractional monomer excitations vanishes 
in this liquid phase, and that its energy spectrum matches that of the Coulomb phase in (3 + 1) dimensional 
quantum electrodynamics. 
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One of the great triumphs of twentieth-century physics was 
to show how the different physical properties of metals, mag- 
nets, semiconductors and superconductors could be under- 
stood in terms of the collective properties of a single ele- 
mentary particle — the electron. For over fifty years, Lan- 
dau's concept of the Fermi liquid, a three-dimensional quan- 
tum liquid whose quasi-particle excitations carry the same 
spin and charge quantum numbers as an electron, has served 
as the "standard model" for metals. Recently, however, 
this theory has been challenged by experiments on a wide 
range of strongly-correlated materials, including quasi-one 
dimensional conductors, cuprate high-temperature supercon- 
ductors and heavy Fermion systems near a quantum critical 
point flU. 

These experiments prompt us to question whether new 
types of quantum liquid, capable of supporting new types 
of excitation might exist in nature? Indeed, excitations 
with fractional charge are known to exist in highly-doped, 
one-dimensional, trans-poly acetylene 01, and in the two- 
dimensional fractional quantum Hall effect 01 . However 
the existence of fractional excitations in three dimensions re- 
mains highly controversial. Indeed, since the understanding 
the quantum Hall effect is bound to the two-dimensional con- 
cept of "anyons" with fractional statistics |5|, it has often been 
argued to be impossible. 

An unambiguous way of resolving this question would be to 
find unbiased evidence of the existence a quantum liquid sup- 
porting deconfined fractional excitations in a concrete, three- 
dimensional microscopic model. This is the goal of this Letter. 
The model we consider is the quantum dimer model 
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FIG. 1: (Color online) Conjectured form for the ground state phase 
diagram of the quantum dimer model on a bipartite lattice in 3D, as 
a function of the ratio p of potential to kinetic energy, following |6, 
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on a three-dimensional diamond lattice, where the first term 
describes the kinetic energy of hard-core dimers tunneling 
between different degenerate configurations on a hexagonal 
"fiippable" plaquette, and fj, sets the ratio of potential to ki- 
netic energy. This model was recently derived as an effective 
description of half-magnetization plateaux in Cr spinels 10,0], 
and also describes spin-less fermions or hard-core bosons on 
a pyrochlore lattice at quarter filling, in the limit of strong 
nearest-neighbour interactions [20]. 

In fact quantum dimer models arise naturally as effective 
models of many different condensed matter systems, and pro- 
vide a concrete realizations of several classes of lattice gauge 
theory § iES 1 ll 11211 . As such, they have become central 
to the theoretical search for new quantum phases and exci- 
tations. A key feature of these models is the existence of a 
gapless "Rokhsar-Kivelson" (RK) point for p, = 1, at which 
all correlation functions exhibit algebraic decay fH). Doping 
the model by removing a dimer introduces two monomers. 
Precisely at the RK point these monomer excitations are in- 
dependent, deconfined excitations, each carrying half of the 
spin/mass/charge associated with a single dimer. 

Field theory arguments suggest that, on bipartite lattices 
in three dimensions, the RK point "grows" into an extended 
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FIG. 2: (Color online) (a) Two adjacent 16-bond cubic unit cells of 
the diamond lattice, showing a dimer configuration in the maximally- 
flippable R-state. The two flippable hexagons contained within this 
picture are shaded blue, (b) Non-local "string" defect (dashed line) 
created by separating two monomer excitations. 



quantum liquid phase — see Fig. Q] Similar behaviour is ex- 
pected in a class of closely related quantum loop models 
3 12]. In both cases the liquid phase is described by the U(l) 
gauge theory corresponding to the Coulomb phase in (3 + 1)- 
dimensional quantum electrodynamics ll, T3, 14, T^. Within 
this effective field theory, monomers are deconfined magnetic 
monopoles, sometimes referred to as spinons, whose mutual 
interactions fall off as 1/r 2 112 111 . Thus, if we can establish 
the existence of a (7(1) liquid phase in the quantum dimer 
model defined by Eq. (Q]i, we will have found a concrete route 
to stabilizing fractional excitations in three dimensions. 

In this Letter we use a range of numerical techniques to de- 
termine the phase diagram of the quantum dimer model on 
a diamond lattice directly from its microscopic Hamiltonian 
Eq. ([T). Our results explicitly confirm the conjectured form of 
the phase diagram shown in Fig. Q] affirming the existence a 
U(l) liquid phase and the absence of a confining string poten- 
tial for fractional monomer (monopole) excitations. 

We begin by considering the ordered phase for fi — ► — oo. 
Here the potential energy dominates, and the ground state 
of Eq. (Q3 is the set of dimer configurations which maxi- 
mize the number of "flippable" hexagons, with one out of 
four hexagons being flippable. This is the so-called R-state 
10], illustrated in Fig. Eta), which has cubic symmetry and 
is eight-fold degenerate. To establish the validity of the pro- 
posed phase diagram Fig. [T] we need to connect this with 
the RK point /j, = 1, for which the ground state the equally 
weighted sum of all possible dimer configurations jjj. We ac- 
complish this using a mixture of exact diagonalization, varia- 
tional Monte Carlo (VMC) and Green's function Monte Carlo 
(GFMC) techniques. 

GFMC is a zero temperature quantum Monte Carlo tech- 
nique which offers a systematic way of improving upon the 
variational wave function output by a variational Monte Carlo 
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FIG. 3: (Color online) (a) Order parameter tur of the R-state as 
a function of /i, calculated using Green's function Monte Carlo 
(GFMC) for [100] clusters with 128, 432, 1024 and 2000 bonds (lines 
and dots serve as a guide for the eye). Results are normalized such 
that 771r = 1 for /j, — > — oo. Exact diagonalization (ED) results 
for 128 bonds are also shown. A sharp jump in rriR can be seen for 
jj,c ~ 0.5-0.7, suggesting a first order transition out of the R-state. 
(b) This coincides with a collapse in the string tension Ai/L. The 
inset in panel (a) shows the finite-size scaling of rriR at the RK point 
{/j, — 1). The inset in panel (b) shows the finite-size scaling of /x c . 
We obtain a value fi c = 0.77 ± 0.02 



calculation |18, 3. Where it converges, GFMC offers accu- 
racy comparable with exact diagonalization 12211 . As an input 
for VMC, we use a trial wave function based on hexagon- 
hexagon correlations, with of order 40 variational parameters. 
Cluster sizes are limited by the rapid growth of the Hilbert 
space with system size — there are 1.3^ dimer coverings of 
an A^-site diamond lattice cluster, and all of these contribute 
to the ground state wave function approaching the RK point. 
However we are able to access families of clusters with edges 
parallel to the [100], [110] and [111] directions which have the 
full (cubic) symmetry of the diamond lattice. Below we fo- 
cus mainly on [100] clusters with 2L 3 diamond lattice bonds, 
where L = {4, 6, 8, 10}. 

In Fig. [3{a) we present numerical results for the order pa- 
rameter ?7ir. This is defined as a sum of projections onto the 
six linearly-independent combinations of the eight degener- 
ate R states. For the smallest 128-bond cluster, wir evolves 
smoothly with ji. There is a perfect numerical agreement be- 
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tween GFMC and exact diagonalization results, and very close 
agreement between these and VMC. For larger systems there 
is a strong suggestion of a first order transition out of the In- 
state — a jump in tor is observed for p c m 0.7 in GFMC 
simulations (2000 bonds). We have studied the finite size scal- 
ing of p c for series of [100], [110] and [111] clusters, and find 
p c = 0.77 ± 0.02 in the thermodynamic limit — see inset 
to Fig. Ob) 1 23 ] . For the clusters considered, tor still takes 
on a finite value for p, > p c , but this is a finite-size effect. 
Exactly at the RK point tor must vanish, and we can use our 
knowledge of the exact ground state to simulate much larger 
clusters. We find that that tor scales to zero as L~ 3 / 2 (see the 
inset to Fig.[3ta)). 

Another indicator of a crystalline, ordered phase is the lin- 
ear, confining "string potential" associated with separating 
two fractional monomer excitations — illustrated in Fig.|2jb). 
Conversely, the energy cost of separating two fractional exci- 
tations by a distance L must vanish in the quantum liquid we 
are seeking. A suitable "string" configuration can be prepared 
by joining the ends of the string defect shown in Fig.HJ?) such 
that it connects opposite faces of a (periodic) cluster of linear 
dimension L. Since GFMC preserves quantum numbers, it 
can be used to calculate the energy of the quantum eigenstate 
corresponding to such a classical string configuration, and the 
"string tension" calculated as Ai/L, where Ai is the energy 
of the string excitation relative to the ground state. 

Deep inside the ordered state, Ai scales as oc pL, how- 
ever the string tension Ai/L collapses abruptly at the value 
of p for which the order parameter jumps — see main panel, 
Fig. (3^b). Taken together, these facts are strongly suggestive 
of a first order transition from a crystaline to a liquid phase 
with deconfined monomer excitations. But they do not yet 
prove the existence of the U(l) liquid we are looking for. In 
order to accomplish this, we must test explicitly the predic- 
tions of the relevant effective field theory. 

The defining property of a dimer model is that every lattice 
site is touched by exactly one dimer. On a bipartite lattice, 
we can implement the this constraint by associating a mag- 
netic field B = VxA with each dimer and each empty bond. 
This field points from site to site along each bond such that 
V • B = at every diamond lattice site B13lll4lll5ll . The total 
flux <f> = J dS ■ B through any plane (cutting bonds) in the lat- 
tice is a conserved quantity. Therefore, for periodic boundary 
conditions, the flux through a set of orthogonal planes defines 
a set of (integer) topological quantum numbers (cf>i, (f>2, <fo). 
We choose to work in a gauge where V ■ A = and in this 
Letter consider only flux sectors of the type (4>, 0, 0), where 
4> = 1 denotes the smallest finite flux and corresponds to the 
state with a single string defect described above. 

This representation clearly has a lot in common with con- 
ventional electromagnetism and, following 11411 . we can use 
this analogy to write down a plausible long-wavelength action 
for the QDM on a diamond lattice 



where E = d t A — V^4o and p± > 0. We have studied p2 
in a perturbation theory about the RK point and find that it 
varies as p2 ~ 0.6(0) x (1 — p). For p > 1, the coefficient 
P2 < 0, and the system chooses those dimer configurations 
with the maximum possible flux <fi — a set of isolated states 
which are not connected by any matrix element of the Hamil- 
tonian Eq. (Q]). All excitations are then gapped. At the RK 
point p2 vanishes, and all flux sectors contribute equally to 
the ground state. The system possesses gapless, transverse 
excitations with dispersion co = p±k 2 . However for p < 1, 
P2 > and the analogy with electromagnetism is complete 
— transverse excitations are now"photons" with dispersion 

W = y/p 2 \k\ Q. 

These photons are the signature feature of the pro- 
posed U(\) liquid state, and offer a beautiful realization of 
Maxwell's laws in a condensed matter system. However, 
Eq. (O also contains information about the finite size scal- 
ing of the ground state energy, as a function of flux (j>. A 
flux <fi through a cluster of volume L 3 corresponds to an av- 
erage magnetic field B = 4>/L 2 . In the £7(1) liquid state, 
this magnetic field is uniformly distributed on the "coarse- 
grained" scale of the effective action Eq. (f2|. It then follows 
from Eq. © that the energy difference A^ = — Eq be- 
tween the ground state of the zero flux sector, and the lowest 
energy state of the sector with flux <fi scales as 
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j d 3 xdt [E 2 - p 2 B 2 - p 4 (V x B) 2 
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Furthermore, by cycling dimers across the boundaries of the 
cluster, we can systematically construct the different flux sec- 
tors <j>. This is the natural generalization of the "thought ex- 
periment" measuring the string tension Ai/L, above. Com- 
bining these results, we have a systematic means of studying 
the spectroscopic signatures of the proposed U(l) liquid state 
in simulations of a finite size systems. 

In Fig. [4] we present this analysis for GFMC simulations 
of an L = 8 cluster of 1024 bonds. Deep within the or- 
dered phase A^ oc L x (f>, and a plot of A^/0 2 at fixed L 
should therefore show a clear distinction between a U(l) liq- 
uid ( A^/(/> 2 oc const.), and an ordered phase ( A^ j<\? oc 1/0). 
For this cluster size a clear division is observed between 
p > 0.7 (liquid) and p < 0.6 (ordered). This is unambigu- 
ous evidence of a phase transition from a linearly confining 
phase into a U{\) liquid. The abrupt, qualitative change in 
the spectra indicates that this phase transition is first-order in 
character. Since we already know the nature of the confining 
phase — the R-state — this completes the phase diagram. 

It is interesting to compare these results with recent work on 
the problem of strongly interacting, hard core bosons at half- 
filling on a pyrochlore lattice lfl7ll . For large values of nearest- 
neighbour interaction, this hard-core bosonic model can be 
approximated by an effective quantum-loop model on the dia- 
mond lattice, closely related to the quantum dimer model stud- 
ied here. Equal-time and static density correlation functions 
for the Bosonic model were calculated using a path integral 
quantum Monte Carlo approach, and an insulating "liquid" 
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FIG. 4: (Color online) Energy gap A^ normalized to the square of 
the "magnetic" flux (j> and distance from RK point, 1 — fi, plotted on 
a log-linear scale for values of fi spanning the ordered R-state and 
the proposed 17(1) liquid phase. In a (7(1) liquid, A^ /<f> 2 ~ const., 
while deep within in the ordered R-state, A^/4> 2 ~ 1/0. Here, a 
clear division is observed between fj, > 0.7 (liquid) and y, < 0.6 
(ordered). Results are from Green's function Monte Carlo (GFMC) 
calculations for a cluster of 1024 bonds. 

state was found, whose finite-temperature correlations were 
well described by the predictions of a U{\) gauge theory. 
However the excitations of this model have yet to be studied, 
and it remains an open question whether this liquid persists 
down to zero temperature. 

In conclusion, the numerical results presented in this Let- 
ter establish the existence of a U(l) liquid ground state in the 
quantum dimer model on a diamond lattice for a finite range 
of parameters 0.7(7) < fi < 1, confirming the field theoret- 
ical scenario proposed by J3, O Q. This liquid state has 
been demonstrated explicitly to support deconfined monomer 
excitations. These results therefore confirm the existence of 
fractional excitations in three dimensions. 

It is our pleasure to acknowledge helpful conversations with 
Kedar Damle, Yong-Baek Kim, Gregoire Misguich, Roderich 
Moessner, and Arno Ralko. We are particularly indebted 
to Federico Becca for advice and encouragement with sim- 
ulation techniques. This work was supported under EPSRC 
grants EP/C539974/1 and EP/G03 1460/1, Hungarian OTKA 
grants K73455 and K62280, U.S. National Science Founda- 
tion I2CAM International Materials Institute Award, Grant 
DMR-0645461, and the guest programs of MPI-PKS Dresden 
and YITP, Kyoto. 



[1] A. J. Schofield; Contemporary Physics 40, 95 (1999). 
[2] C. M. Varma, Z. Nussinov, W. v. Saarlos, Phys. Rep. 361, 267 
(2002). 

[3] W. P. Su and J. R. Schrieffer, Phys. Rev. Lett. 46, 738 (1981). 

[4] R. B. Laughlin, Phys. Rev. Lett. 50, 1395 (1983). 

[5] F. Wilczek, Phys. Rev. Lett. 49, 957 (1982). 

[6] D. L. Bergman, et al, Phys. Rev. Lett 96 097207 (2006); ibid 

97, 139906 (2006) 
[7] D. L. Bergman, G. A. Fiete, and L. Balents, Phys. Rev. B 73 

134402 (2006). 

[8] D. S. Rokhsar and S. A. Kivelson, Phys. Rev. Lett. 61,2376 
(1988). 

[9] E. Fradkin, Field Theories of Condensed Matter Systems, Ad- 
dison Wesley, New York (1991). 

[10] R. Moessner and S. L. Sondhi, Phys. Rev. Lett. 86, 1881 (2001). 

[11] G. Misguish and C. Lhuillier, in Frustrated Spin Systems, Ed 
H. T. Diep, World Scientific, Singapore (2004). 

[12] F. Trousselet, D. Poilblanc and R. Moessner, Phys. Rev. B 78, 
195101 (2008) 

[13] D. A. Huse et al, Phys. Rev. Lett. 91, 167004 (2003). 
[14] R. Moessner and S. L. Sondhi, Phys. Rev. B 68,184512 (2003). 
[15] M. Hermele, M. P. A. Fisher, and L Balents, Phys. Rev. B 69, 
064404 (2004). 

[16] P. Fulde, K. Penc and N. Shannon, Ann. Phys. (Leipzig) 11, 892 
(2002) 

[17] A. Banerjee et al. Phys. Rev. Lett. 100, 047208 (2008). 
[18] N. Trivedi and D. M. Ceperley, Phys. Rev. B 41, 4552 (1990). 
[19] M. Calandra Buonaura and S. Sorella, Phys. Rev. B 57,11446 
(1998). 

[20] In this special case, the Fermi sign can be "gauged" out of the 
problem by a suitable choice of convention for labeling lattice 
sites. 

[21] For classical analogues to this problem, see e.g. C. Castelnovo 
et al, Nature 451 42, (2008); G. Misguich et al. Phys. Rev. B 
78, 100402(R) (2008) 

[22] We have checked our GFMC results explicitly against exact di- 
agonalization for small system sizes, and for large systems sizes 
against an expansion in 1/|/^| for /i — > — oo, and perturbation 
theory in 1 — /i for fj, — > 1. We find complete agreement in all 
cases. This (lengthly) analysis will be presented elsewhere. 

[23] The L~ 4 -scaling follows from the fact that the R-state is 
gapped, while the competing liquid phase has linearly dispers- 
ing photon excitations. 

[24] We work in units such that the volume of the 16-bond cubic unit 
cell is equal to 8. This, together with the normalization chosen 
for the the flux <j>, defines the units of p2 and therefore the speed 
of light for photon excitations. 



